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The emergence of gauge particles (e.g., photons and gravitons) as Goldstone bosons arising from 
spontaneous symmetry breaking is an interesting hypothesis which would provide a dynamical set¬ 
ting for the gauge principle. We investigate this proposal in the framework of a general SO (N) 
non-Abelian Nambu model (NANM), effectively providing spontaneous Lorentz symmetry breaking 
in terms of the corresponding Goldstone bosons. Using a nonperturbative Hamiltonian analysis, 
we prove that the SO (N) Yang-Mills theory is equivalent to the corresponding NANM, after both 
current conservation and the Gauss laws are imposed as initial conditions for the latter. This equiv¬ 
alence is independent of any gauge fixing in the YM theory. A substantial conceptual and practical 
improvement in the analysis arises by choosing a particular parametrization that solves the non¬ 
linear constraint defining the NANM. This choice allows us to show that the relation between the 
NANM canonical variables and the corresponding ones of the YM theory, A“ and A 6 -' , is given by 
a canonical transformation. In terms of the latter variables, the NANM Hamiltonian has the same 
form as the YM Hamiltonian, except that the Gauss laws do not arise as first-class constraints. 

The dynamics of the NANM further guarantees that it is sufficient to impose them only as initial 
conditions, in order to recover the full equivalence. It is interesting to observe that this particular 
parametrization exhibits the NANM as a regular theory, thus providing a substantial simplification 
in the calculations. 
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I. INTRODUCTION 

The possible interpretation of gauge particles (e.g., photons and gravitons) as Goldstone bosons (GBs) arising from 
some spontaneous symmetry breaking, dates back to the pioneering works ofi and HI- The former used the 
standard coset construction of the effective theory jdj, in which the spontaneous Lorentz symmetry breaking (SLSB) 
is realized nonlinearly in terms of the GBs, with the matter fields transforming linearly under the unbroken subgroup. 
Their conclusion was that any gauge theory is a theory of some spontaneously broken symmetry. A recent application 
of the coset construction to general relativity was reported in Ref. jsj. 

On the other hand, Refs. [2, '3j considered the explicit case of spontaneous Lorentz symmetry breaking to produce 
the photon field as the condensate arising from a self-coupled four-fermion model, following similar steps as those 
developed in Ref. Q to describe the superconductor solutions in field theory. Similar ideas have been revisited for 
photons and also have been extended to gravitons in Refs. |7H9|| . 

An alternative approach was proposed by Nambu in Ref. |10j . where the emphasis was shifted to the description of 
the SLSB system (QED, in this case) only in terms of the GB degrees of freedom (d.o.f.), which were introduced via a 
nonlinear constraint, similarly to the nonlinear sigma model description of pion interactions. Such a cr-QED is defined 
by the Maxwell’s Lagrangian plus the constraint A^A^ = n 2 M 2 , which is to be substituted into the Lagrangian. 
Here, is a properly oriented constant Lorentz vector, while M is the proposed scale associated with the SLSB. 
This constraint can be understood as providing a nonzero vacuum expectation value (A^) = n^M. Nevertheless, the 
goal in this model was to show that it is in fact equivalent to standard QED, instead of yielding a physical violation 
of the Lorentz symmetry. This equivalence was manifest up to the tree-level calculations studied in Ref. 0- Later 
on, these calculations were extended to some processes at the one-loop level, with identical results: all contributions 
arising from the SLSB sector of the model canceled out, yielding the standard QED results 0- The cr-QED model 
has been further studied [HI and extended to the non-Abelian [l3l - [l6| and gravitational cases 0 , which we will 
generically call generalized Nambu models. General conditions regarding how the gauge symmetries were recovered 
from the corresponding SLSB models are worked out in [ij§-[20|. 

Perturbative calculations in the non-Abelian case show again that, to the order considered, all SLSB contributions 
to physical processes cancel out, yielding an equivalence with the starting Yang-Mills (YM) theory, in complete 
analogy with the Abelian case. This fact has been interpreted by stating that the corresponding nonlinear constraint, 
which defines each Nambu model, can be interpreted as just a gauge choice in the associated Abelian or non-Abelian 
gauge theory. This would lead to an equivalence between the Nambu model and the corresponding gauge theory in 
a fixed gauge. Nevertheless, this statement requires some qualifications, (i) To begin with, the number of d.o.f. of 
the Nambu model is larger than that of the corresponding gauge theory, which can be understood because the former 
has lost gauge invariance, (ii) Fixing the gauge in any gauge theory requires the introduction of ghost particles (via 
the Becchi-Rouet-Stora-Tyutin procedure, for example Ref. 0 ) which play a fundamental role as internal particles 
in calculating physical processes. Thus, in order to establish the proposed equivalence, one would need to study the 
contributions of the ghosts to physical processes. A possible decoupling of them is by no means evident, especially 
due to the nonlinear character of the proposed gauge fixing. The first point (i) has been taken into account in most 
previous works and was emphasized in Ref. [22:]. The general statement, phrased in different ways in different papers, 
is that the Nambu model is equivalent to the corresponding gauge theory only after current conservation together 
with the Gauss laws have been imposed as initial conditions, since the dynamics of the Nambu model preserves their 
conservation for all times. The second point (ii) has not been considered at all and will be dealt with in a separate 
publication, for the case of the Abelian Nambu model [23|. 

Recently, the study of possible observable violations of Lorentz invariance has attracted considerable attention, from 
both the experimental and theoretical points of view. Explicit Lorentz symmetry violation is found to be incompatible 
with the Bianchi identities fH | and therefore this approach is not consistent with general relativity, unlike SLSB, 
where this issue does not occur. The construction of the Standard Model extension performed by Kosletecky and 
collaborators [Ht is a framework in which Lorentz violation is considered as arising from a spontaneous symmetry 
breaking in a more fundamental theory. A distinguished class of models in that framework are the so-called bumblebee 
models, which are tensor theories exhibiting physical SLSB. They include GB modes, and depending on the explicit 
form of the theory they have additional modes and constraints. These models have been thoroughly investigated 
in relation to electrodynamics [2^, 1H3] and gravity [28j - [3l| . As a matter of fact, generalized Nambu models can be 
thought of as very particular cases of bumblebee models, where the non-Goldstonic d.o.f. of the latter are frozen, 
leaving only the GB excitations. 

In the present paper we generalize to the non-Abelian case and improve the nonperturbative Hamiltonian analysis 
developed for the Abelian Nambu model (ANM) in Ref. [22[. The non-Abelian Nambu model (NANM) associated to 
the group SO(iV) is defined by the YM Lagrangian 

£(A“) = -\f; v F^ - 


(1) 
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plus the condition 


A“A“ M = n 2 M 2 , M 2 > 0, /x = 0,1,2,3, a = l,...,N, (2) 

which is to be solved and substituted into the YM Lagrangian. Our goal is to determine which additional conditions 
have to be imposed upon the NANM in order that its Hamiltonian reduces to that of the YM theory. No discussion 
is provided of the possible perturbative equivalence of the so-corrected Nambu model and the YM theory in the fixed 
nonlinear gauge. 

We proceed via the following steps, (i) We start by constructing the Hamiltonian for each type of NANM (depending 
on n 2 it can be time-like, space-like, or light-like) in terms of the corresponding canonical variables for each case, (ii) 
We show that these Hamiltonians are related via a canonical transformation to a Hamiltonian that has the same form 
as the YM Hamiltonian in the standard variables Af, Ef , i = 1,2, 3, except for the fact that the Gauss laws = 0 
do not appear as constraints; nevertheless, the canonical transformation leads to the correct brackets between the 
canonical variables Af,Ef arising from the canonical algebra of each NANM. (iii) We prove that the NANM dynamics 
preserves the evolution of 0 ° in such a way that it guarantees that the imposition of fi° = 0 for some initial time leads 
to Q a (t) = 0 for all times. In this way, enforcing the Gauss laws as first-class Hamiltonian constraints at some initial 
time makes the NANM equivalent to the corresponding YM theory in a nonperturbative way and independently of 
any gauge fixing. Consistency with the NANM dynamics avoids the generation of additional constraints. 

The present paper is organized as follows. In Sec. II we consider the specific case of the space-like NANM (SL- 
NANM) by solving the nonlinear constraint (© in terms of Affz 3 and starting with the remaining 47V— 1 d.o.f. per point 
in coordinate space. The Lagrangian equations of motion are obtained and the canonical momenta together with the 
canonical Hamiltonian are subsequently constructed. The standard variables of the YM theory Af and E b i ,i = 1,2,3 
are written in terms of the canonical variables of the SL-NANM and their induced algebra is calculated, which is 
summarized in Appendix B. The analysis of the SL-NANM in terms of the Dirac method reveals that this model has 
additional second-class constraints, which are further imposed strongly by introducing Dirac brackets to obtain the 
d.o.f. of the reduced phase space, together with their algebra. The induced Dirac-brackets algebra for the variables 
Af and E b i are also calculated. The final extended Hamiltonian for the SL-NANM, rewritten in terms of the variables 
Af and E b i, is finally obtained and compared with the standard YM Hamiltonian. The conditions under which both 
theories are equivalent are established, which requires understanding the time evolution of the Gauss functions Q a 
under the SL-NANM dynamics. A similar analysis can be carried out for each of the remaining cases, corresponding 
to the time-like and light-like NANM. Section III presents a substantial conceptual and practical improvement over 
the previous individual calculations to study the relation between the NANM and YM theories. We start from a 
general parametrization that solves the constraint © for arbitrary values of n 2 , in terms of 3 TV d.o.f. $( 4 , A = 1,2, 3, 
and repeat the canonical analysis, paying attention to the relations between the standard canonical variables of the 
YM theory (Af, E hj ) and those of the NANM (<£( 4 , n^). The calculations are enormously simplified after one realizes 
that the transformation ($( 4 ,^) —»■ (Af,E b ^) is a canonical transformation, once the E b i are recognized as the 
momenta canonically conjugate to the Af via the kinetic part of the NANM Hamiltonian action. Another useful 
property of the new parametrization is that it exhibits the NANM as a regular theory (i.e., no constraints appear in 
the Hamiltonian analysis). This is proved in Appendix C. The canonical NANM Hamiltonian, rewritten in terms of 
the YM variables (Af,E b i), is finally obtained and the conditions under which it reduces to the YM Hamiltonian are 
determined. Again, this requires the calculation of the time evolution of the Gauss functions 0“ under the NANM 
dynamics. Finally, we close with a summary and some comments in Sec. IV. Appendix A serves to establish notation 
and briefly reviews the canonical version of the SO (TV) YM theory, which we use as a benchmark to identify the 
conditions under which it is equivalent to the different realizations of the NANM. 


II. THE SPACE-LIKE CASE OF THE NON-ABELIAN NAMBU MODEL (SL-NANM) 

Before considering the specific case of the SL-NANM, let us recall some general properties of the NANM. The 
Lagrangian is 

C N ANM(Af) = + A (AfA a » - n 2 M 2 ) , M 2 > 0, a = 1,..., TV, (3) 

with the notation and conventions introduced in Appendix A. Here, A is a Lagrange multiplier and the vector is 
such that n 2 = ±1,0. 

The general procedure by which we will analyze each model is to explicitly solve the condition 

Af,A afi = n 2 M 2 

A" 


(4) 
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and substitute the adequate parametrizations directly into the Lagrangian £nanm(^), defining in this way the 
canonical degrees of freedom of the model. Subsequently, we obtain the corresponding Hamiltonian and determine 
its relation to the YM Hamiltonian (IA28I) together with the canonical algebra (IA29I) . As expected, the equations of 
motion of the NANM will not be those of the Yang-Mills theory arising from £( A M ) in Eq. (IA1I) . This property will be 
explicitly shown in the remaining sections of the paper. In this way, the conservation of the current J“ does not follow 
as a consistency condition from the equations of motion in the NANM, as happens in the YM case. We will show that 
the canonical structure of the NANM will induce the standard-algebra YM [ (IA29I) ] together with a Hamiltonian that 
differs from Eq. (IA28I) by the property that the Gauss laws do not appear as constraints. Nevertheless, the dynamics 
of the NANM guarantees their validity for all time, once they are imposed as initial conditions. 

The standard solutions of the condition Ag A a/i = n 2 M 2 , arising from the different choices of n 2 , are 


> 0 : Ai = JM 2 + - A b 0 A b 0 , 6 = 2,3,.., TV, 


n 2 < 0 : A\ = \JM 2 + (Ag ) 2 - (A “) 2 - (Ag ) 2 - (Ag ) 2 

A\A b \ 

4 B 2 ) ' 


n 2 = 0 : Ag = B a ( 1 + 


^ = B “ ( 1 - T0) 


* = 1 , 2 , 


(5) 

( 6 ) 
(7) 


which define the NANM in its time-like (TL-NANM), space-like (SL-NANM) and light-like (LL-NANM) represen¬ 
tations. In the time-like and space-like cases we start with 47V — 1 d.o.f. per point, while in the light-like case this 
number is 3TV. Next we concentrate on the SL-NANM. 


A. The equations of motion in the SL-NANM 


We start with the extension to the non-Abelian case of the parametrization A 3 = y/M 2 + Ag — A 2 — A 2 , which 
is frequently used in the Abelian case to exhibit the remaining SO(2,1) symmetry after the spontaneous Lorentz 
symmetry breaking. The Lagrangian constraint (jTj is now solved for Aglig. In this way we start from 47V — 1 d.o.f. 
in coordinate space, which we denote in the following way: Ag, A|, Aj5, with a = 1,2,...,TV; a = 2,3, k = 

1, 2,3; k = 1, 2. The numbers of each of the corresponding fields are TV, 2, 3(TV — 1), respectively. With this notation, 
we write 


AJ = sjM 2 + (Ag ) 2 - (Al ) 2 - (At) 2 , 


( 8 ) 


which exhibits the remaining symmetry group SO (TV, 37V — 1). As a matter of notation, superscript indices label a 
group index, while subscript indices refers to a space-time index. 

In the notation of Appendix A, the Lagrangian density ([3]) takes the form 


£(Ag, AlA%) = l ( (Ef ) 2 + (El ) 2 + (Eg ) 2 - B%B%) - J 0 “A“° + J\A\ + J?Af + j\A\. 


(9) 


The equations of motion are 


A a 

S ia_ g 31i =0) 

A 3 


( 10 ) 


with the notation 


rOa 


A a 

+ £• 31110-0 
+ Al" 0 ' 


s va = - r) a . 


(ii) 


( 12 ) 


The numbers of equations in Eq. m is only (3iV — 1) because the simultaneous choice 2 = 3 and a = 1 does not 
appear, as is a function of the dynamical variables. 

Since Ag is just a shorthand for ©, it follows that 


AS . A- . A a . 

\1 ^0 \a _ k \1 \a 

3 “ (4 1 0 A - s 7 1 ‘ 

3 ^3 Ai 


(13) 
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The fields E\ and E a are given by 


E\ =Al- D 3 Al = 


A a A- A a \ 

A a _ k A 1 _ _L A a I _ Do A 1 

Al A Q 4 l /1 fc A lAi I U 3 A 0i 
k ^3 ^3 ^3 / 

Et = M - DiAl, 


E\=A\-D- k Al 

which, for the moment, constitute a compact way of identifying the velocities Aft , Af, and A\. 


(14) 


B. The Hamiltonian density of the SL-NANM 

The canonically conjugate momenta are 


TT a — 
Ai o — 

dC 

oAq 

A a 

— 

(15) 


dC 

. ,A\ 


n lfc = 

dA\ 

~ h k h 3-^ T» 

(16) 


dr 

, A a 

(17) 

n ai = 

d if 

_ TACL 771 1 l 

— -^3 ^1 5 


with nonzero Poisson brackets (PBs) 

{Ag(x,t),ng(y,t)}=<5 3 (x-y)r b , {4(x,t),II lS (y,t)} = <^ 3 (x-y), {4(x,t),I#'(y, t)} = <^ 3 (x-y). 

(18) 

Next we solve for the velocities. From Eqs. dm), (hhd and we (respectively) obtain 


A a A- A a \ ttI 

_ib _ ilk A- — —A a I _ Dr, A 1 — — 4 1 

Al A 0 41 A k Al A l a1 A 3 , 

k A 3 a 3 a 3 ) ^0 

PI 1 

A a — n A a - llfi A a — TT al 

A i U i A 0 Al A i ~ 11 ■ 

^0 


From (BUI) and (BT1) we can solve for 


IP 


Ai = tf k + D- k Al + -EAl 

A o 


n 1 

il 0 A a 


A a = U a> + D . A a + U 
■A-n 


(19) 

( 20 ) 
( 21 ) 

( 22 ) 

(23) 


We can substitute these velocities into Eq. dm but we cannot solve for all the Ag which enter into the sum AgAg. 
At most we could solve for one velocity, say Ag, in terms of the remaining (TV — 1) Ag. This is consistent with the 
existence of (TV — 1) primary constraints 4>g, which we choose as 


n 1 

d> a — TT a _0 A a 

V 1 “ li 0 41 ^0 ’ 
^0 


(24) 


arising from Eq. m It is more convenient to consider the solved velocity as A 3 , encoded in the definition of E 3 and 
written in terms of the canonical variables as 


i Ai 

ri 3 


£3 =IP —, 
^0 


(25) 


via the remaining relation m corresponding to a = 1. 

From Eqs. (USD, (USD and CH> we can express the electric fields in terms of the canonical momenta as 


E 1 - A 1 l — 

^3 ~ a 3\ 


E\ = n H + 


e° = n Q 


55 ' U a 
Ab) Ai ’ 


(26) 
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where we have used Eq. (l25l) . 

In Appendix B we show that the above definitions of the electric fields Ef in terms of the canonical momenta 11“, 
together with the canonical algebra (fTSl) lead to the following PB relations 


A b j(y,t)} = 0, {E al (x,t),E b:1 (y,t)} = 0, {A“(x, t), E b: > (y, t)} = ~6jS ab S(x - y), (27) 

which reproduces the YM algebra (IA29I) . In the following we will also need the PB of the variables Ag, ng with A“ 
and Ef, which are summarized in Eqs. (1B9I) and (IB 101) of Appendix B. 

The next step is to obtain the Hamiltonian density 

n = ngAg + n H ii + n ak Af - L(Ag, a?), ( 28) 


where C(Af,A\,Af) is given in Eq. ([9]). The main goal is to rewrite this Hamiltonian in terms of the fields Ef and 
Aj in order to compare with the YM result (IA28I) . From Eqs. (1151) . (flTH) and (fT71) . we substitute the momenta into 
the above equation obtaining 


n = El^Al + (££ -E^j Ai + - E^EE j A% - L(Ag, Af), 

A? 


A} 


A' 

pl 


= E\ 


A a . A 1 

jjO A a _ ELaI _ 

ai a 0 a 1 A % 

^3 ^3 


_ k Aa 

A\ Ak 


a 3 , 

7i a Aak 


+ EfA H + EfA ak - L(Ag, A?) 


(29) 

(30) 


where we subsequently replace the velocities Af and Ag in terms of the electric fields, using Eqs. m ■ This yields 
E = E\ {E\ + D- k Al ) + Ef {Ef + D z Af) + E\ {E\ + D 3 Ag) 


-- ( (Ef) 2 + (Ef) 2 + (E\) 2 - BfBf) + JqA a0 - JlA\ - jf Af - J\A\ 


(31) 


which can be further rearranged as 


1 


1 


1 


U = -EfEf + -EfEf + -EfEl + -BfBf + JfA M - JfAf - JfA “ - JfA 


k ^k 


+E\D k Al + EfD t Al + ElD 3 Al, 
n = X -EfEf + l -BfBf - JfAf - Ag (D z Ef - J °). 


(32) 


Here we have integrated by parts the term Ef (DiAf) in the Hamiltonian. The form of the above Hamiltonian density, 
together with the PBs in Eq. (1271) . is similar to that of the SO (IV) YM theory (IA28I) . except for the following facts: 
(i) the coordinates Ag are dynamical instead of being Lagrange multipliers, (ii) the Gauss functions, defined as 


n a = (DiEf - J 0 a ) = £' 


0 a 


(33) 


are not constraints in the SL-NANM and (iii) we have the additional primary constraints E71) . 

Then we need to continue the Hamiltonian analysis by applying Dirac’s procedure starting from the extended 
Hamiltonian density 


He = \EfEf + l -BfB a k - JfAf - AgD“ + f 


n 1 

TT“ _ A a 

n o cA 


The evolution of the primary constraints yields 




J d 3 y {^E b E b + X -B a k Bl - JfAf - Ag (D,Ef - J“)) J , 
|d>“, J d 3 y (—Aj^-Agn" 6 )}, 


the calculation of which requires the following PBs calculated in Appendix B 

{$?, Ag}=0, {3>lE b }= 0, {$;, i?n=o, 

{*!> Aj } = ^|, {<&?, Ag } = -S si . 


(34) 


(35) 


( 36 ) 
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In this way, the only contribution arises from the terms proportional to Ag in the Hamiltonian density. The result 

4' 


$? = - + s ab n b 


produces secondary constraints, which we write as 


a 1 

<f) a — A a _ 0 pg 

®2-A 0 q1 U . 


Next we calculate the time evolution of $2 using the relations 

{fi s , <I>f} = 0, {$“, Al} = -6 5 

which are included in Appendix B. We obtain 


rah 




n a Ag 
fYAg’ 


(37) 


(38) 


(39) 


= W a + n a + 


A a A 1 

Ain AT. 


0^0 e 
2 " ■ 


K) 


In fact, we can solve 


404 c 

<5 as + 

v K) 2 , 


^ = -w s 


for the arbitrary functions /r°, concluding that 


l l a = _ I S ac_ 


(a^ 0 ) asac °J wc - 


In this way the Dirac method stops and we are left with 2(N — 1) constraints 

n 1 - - a 1 

(h a — TT a llO Aa <T>a _ Aa ^0 pa 

— lio A lA °’ ^2—^0 Qi Si ) 

which are second class. Thus the number of d.o.f. per point of the SL-NANM is 

#d.o.f. = 1(2(4 N - 1) - 2 (N - 1)) = 3 N, 


(40) 


(41) 


(42) 


(43) 


which does not correspond to the number of d.o.f. of the SO(N) Yang-Mills theory. 

The next step is to set the constraints C^ll strongly equal to zero, in order to eliminate the variables Ag and ng, and 
to subsequently introduce the corresponding Dirac brackets among the remaining variables. To this end we require 
the matrix constructed with the PB of the constraints 


M = 


R ah 

rpba 


rj^ab 

Sab 


R T 
-T t S 


where 


R 3 ^ = |<f> a , T aa = {d> a , $ a }, S ,a5 = {$g, $ a } . 

The required calculations produce 

{$?, $ 1 } 

'Al 


(44) 


(45) 


R ab = Q) T ab = 


Ag Ag 


= - \S ab + 1 = T 


’-iba 


A J 


S ab = 


0 

D 1 


' - O b - O a 

s~iabm | /~ilarri s-ilbm 

6 +c m 


(AA") ■ 


(46) 
















according to the results in Appendix B. 
The matrix T is invertible, yielding 


(T-i y ab 



KK ) ’ 


(47) 


in such a way that 


M -1 


T~ 1 ST~ 1 

T-i 


-T-i 

0 


The Dirac bracket is 

{A(x),B(y)}* = {A(x),B(y)} - {A,4>\}(T~ l ST^f 1 B} 

+{A, mT-Y b {4 L B} - {A, ^IKT- 1 ) 55 ^, B}, 

which leads to the result 


(48) 


(49) 


{A(x), B{y)}* = {A(x), B(y)}, (50) 

for the YM variables A“ and E hj . The above conclusion arises from the fact that each of the additional PBs in Eq. 
(l49l) includes a contribution from which has zero PB with those variables, according to Eq. (IB 111) . In other words, 
we recover the algebra 

{A?(x,i),A$(y,i)}* = 0 = {£“(x,t),£«(y,t)}*, {A“(x,f),£^(y,t)}* = -S{5 ab 6(x - y), (51) 

corresponding to the YM theory given in Eq. (IA'291) of Appendix A. Having set the constraints (1T21) strongly equal to 
zero, the extended Hamiltonian (1M1) now reduces to 


U E = \e?E? + l -B a k B a k - J?A“ - Agfl°. 


but we are still missing the Gauss laws O a = 0, because the Ag are dynamical degrees of freedom. 


(52) 


C. The evolution of the Gauss functions f l a in the SL-NANM 

Next we study the time evolution of the Gauss functions, starting from the Hamiltonian density (ld2l) . A direct use 
of (IA29D leads to 


= - gC abc A b o n c - D^ a - J d 3 y {fi a (z), A b 0 (y)} n b (y). (53) 

where f l a (x) —> (6 ac dk x + gC abc A b k ) E k inside a PB because J 0a has been considered as an external current. Since 
{Ag, A k } =0, we only need the brackets {E k (- r), Aq(j/)}. Using Eq. (IB9I) we obtain 

n a = - g C abc A b n c - D k J^ a + D k . (54) 

The above equations guarantee that by (i) imposing current conservation D k J lia = 0 at some initial time t = 0 and 
(ii) demanding that the Gauss laws fU = 0 hold at t = 0, we obtain d 0 Q a = 0 (a = 1,2,..., N) as well at t = 0. 
This is enough to prove that with these two initial conditions, the Gauss laws will hold for all time. In this way 
we can recover the SO (N) Yang-Mills theory by imposing the Gauss laws as Hamiltonian constraints, with arbitrary 
functions N a adding —N a Q a to He and redefining Ag + N a = O a . This leads to 

He = ^(E 2 + B 2 ) - O a tt a + J t a A ia , (55) 

where the O a are now arbitrary functions, and thus we get back to the YM Hamiltonian density (|A28j) . The subsequent 
emergence of the SO (At) YM theory guarantees current conservation for all times, as a consequence of the equations 
of motion. 
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From the perspective of the GB modes, the situation in the SL-NANM is as follows. We have started from a theory 
invariant under SO (TV, 37V) defined by Eqs. 0 and (0. Solving the constraint 0 in terms of A 3 means that we 
have the nonzero vacuum expectation value (A 3 ) = M , which breaks the symmetry down to SO(7V, 3TV — 1), with the 
appearance of (47V(47V — l)/2 — (4iV — 1)(47V — 2)/2) = 47V — 1 GBs. Nevertheless, the SL-NANM phase space still 
contains 2(7V — 1) second-class constraints [Eq. (l42l) ] which can be imposed strongly, yielding a reduced phase space 
with i(2(4iV — 1) — 2(TV — 1) = 37V coordinates per point. A similar analysis yields 3TV as the number of independent 
GB modes in each realization of the NANM. This conclusion is consistent with Ref. [l5|, where the final independent 
GB modes are denoted by a 1 , (i = 1 // = 1,2,3) for the time-like case and a 1 ,, (i = 1, ..,1V, \j!' = 0,1,2) 

for the space-like case. Once more, we verify that in order to regain the 27V independent vector bosons of the YM 
theory, we still have to impose the TV Gauss laws Q a = 0 as first-class constraints. In this way we are left with 
2 (61V — 21V) = 2TV coordinates per point. 

To summarize, the emergence of the SO(IV) YM theory from the SL-NANM can be established only after imposing 
both current conservation and the Gauss laws as initial conditions. 


III. A UNIFIED DESCRIPTION OF THE NON-ABELIAN NAMBU MODELS 


A procedure similar to that presented in the previous section for the SL-NANM can be repeated for the TL-NANM 
and the LL-NANM, with identical results. The standard variables A“, E b of the YM theory can be expressed in terms 
of the canonical variables of each version of the NANM, the algebra of which induces their brackets to be those of Eqs 
(IA29I) . These transformations also allow the canonical Hamiltonian density of the NANM to be rewritten in terms 
of Af and E b , yielding a result with exactly the same form as Eq. (IA28I) . Nevertheless, there is a main difference 
between the so-constructed Hamiltonian and the full YM Hamiltonian. In fact, the Gauss laws fl a = 0 do not appear 
as constraints in the NANM (since the Ag are not Lagrange multipliers), but rather as functions of the respective 
d.o.f.. 

This has motivated us to search for a unified and simpler discussion of the generic NANM. To this end, we find it 
convenient to generalize the parametrization 0 to all cases in the form 

A o=B a (l + ^), A a 3 = B a (l-^J, N = (A b A b + n 2 M 2 ), AB 2 ± N ? 0, (56) 


which certainly satisfies the condition 0 and is written in terms of the 37V independent GBs ( B a and A|). 

The relation between the above parametrization and the purely Goldstonic d.o.f. introduced in Ref. [151 ]. which 
relabel as a b , can be established as follows. The 47V fields a b are subjected to the TV additional constraints 


we 


= 0 , 


(57) 


leaving only 3N independent GB modes. In terms of them, the original fields are written as 

■>b 


n,, 


K = a ^~l ( M2 “ n 2 a 2 ) 1/2 , n£ = n M s\ s 2 = 1 , n 2 ± 0 , 


(58) 


which satisfy the condition 0. Equation (1551) can be inverted to produce 




(59) 


which allows us to express a b = a b {B c , A“) by employing Eq. (|56)l . 


A. The equations of motion 

After the substitution of (1561) into the Lagrangian density 0, the variation of the corresponding action with respect 
to Ag yields 


0 = J d i x(D /1 F 111 ' - J u ) a 5A a v , 


( 60 ) 
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where the 5A “ are not all independent. In our case, Eq. m leads to 


^o“ = 


! ,N_\ tab N 2 B“B b 
4B 2 J 4B 2 B 2 


SB b +^A*6Al 


SA a 3 = 


-JL\ 


i - — i s ab 
4B 2 J 


N 2 B a B 


b 1 


SB” ~ 


4 B 2 B 2 

in terms of the independent variations SA b and SB a . In this way the equations of motion are 

r>b 

6A“ : + [£° b - £ 3 1 A " = 0, 


(61) 


(62) 


(63) 


6B a : 



N 2 B b B° \ ob 
4 B 2 B 2 ) 



N 2 B b B* \ Zb 
4 B 2 B 2 ) ’ 


(64) 


in the notation of Eq. (USD- 

Let us recall that in the case of the SO(IV) YM theory the equations of motion are just given by £ va = 0. Also, the 
above equations of motion do not imply current conservation D„J va = 0 , basically because the condition 0 breaks 
non-Abelian gauge invariance. A way to recover the YM equations of motion together with gauge invariance is to 
impose the Gauss laws £ 0a = 0. In this way, under the conditions 4 B 2 ±N Eq. (1M1) yields the solution £ 3b = 0. 
These two conditions in Eq. (1631) provide the final set £ la = 0. 


B. The Hamiltonian density 

In order to unify the notation when going to the Hamiltonian formulation we introduce the 3AT d.o.f. 4>^, A = 1,2, 3, 

= = (65) 

in such a way that the coordinate transformation 

A i = ( 66 ) 

arising from Eq. (15611 is invertible. In fact, the inverses are 

AS 


<!>" = Aj\ $“ = A“ 


$3 = 


2\/Apf 


\f A ^3 + \f A W 


+ n 2 M 2 . 


We also have 


A“ = Ag(4>$t), 


(67) 


( 68 ) 


in terms of Eq. (l 66 l) . according to the first relation in Eq. (l56l) . The relevant property of the transformation (l 66 l) is 
that 


f) A a 

M = -Tzr?-^ 


19A“ 9A“ 


84> b B 


d<S> b B 84> b B ’ 


together with the invertibility of the velocities 


cM>“ .. 

= A A b 
A dA b *' 


(69) 


(70) 


In the following we will not require the explicit form of the transformations (1561) and (1671) . but only their generic 
form (1661) . together with the property that this transformation can be inverted. 

Next we proceed to calculate the Hamiltonian density of the NANM in terms of the canonically conjugated variables 
$(4 and n^, and we employ a procedure that allows us to make direct contact with both the YM Hamiltonian density 
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(IA28I) and the YM canonical algebra (IA29I) . After making the substitutions (l66l) and (l68l) . the Lagrangian density 
(131) can be rewritten as 


where 


4anm($,$) = - \b, ?b? - J a >‘A; 


Et=A^D t A a 0l B? = -e ijk F“ k , 


(71) 


(72) 


with E? = _£“(<!>, 4) and B"' = H“(<f>). The canonically conjugate momenta are calculated as 

_ 5 £ nanm (^; $) = F b 9A b = b dA\ 

A 54“ *54“ ‘5$^’ 

employing (l69l) . The inverse of Eq. (l66l) allows us to write the colored electric fields Ef as functions of the momenta 
of the NANM 


The Wronskian of the system is 

5 2 £nanm($, 4) 


det 


54“54^ 


Ei(*, n) 


d<S> a A 


n 


a 

A • 



= det 


( dA? 5aA 

\d* b B d$ a A ) 


= det 


( dA? 5A“ \ 

\d$ b B d$ a A ) 


*0, 


(73) 


(74) 


as shown in Appendix C. In this way, the NANM is exhibited as a regular system in the parametrization (1531) . so that 
no constraints are present. 

The NANM Hamiltonian density is 


Hnanm = n a A®A - (\e°E« _ - J^aA , (75) 

which we rewrite in successive steps 

f)(b a /1 i \ 

Hnanm = n a-q^A - l^ E ?E? - 2 B ?B? - J 0 M£j , (76) 

Hnanm = E\A\ - Qe“E“ - , (77) 

Hnanm = E\ (E b + DiA b 0 ) - Qe“E“ - , (78) 

H N anm( 4>, n) = Ef + ±B“B“ - (DiE\ - J b0 ) A* + J“A“, (79) 

where we have used Eqs. GOD, USD, and (1731) . together with an integration by parts in the term containing the 
covariant derivative. The dependence of Hnanm on the canonical variables $,n is clearly established by the change 
of variables (BSD, (BSD and (USD- The canonical variables of NANM satisfy the standard PBs 

{<^(x),d> fc B (y)}=0, {H a A {x),H b B {y)} = 0, {^(x),H^(y)} =6 ab 6 AB 6 3 (x-y). (80) 


Now we can consider the NANM Hamiltonian density (1731) from the perspective of the fields A“ and H“. The relation 
arising from the velocity-dependent term of the NANM Hamiltonian action, 


d A x W A $ a A = / d A x E?A“ = / d A x (-E“)A“ 


(81) 


[used previously in obtaining Eq. (1771) ]. establishes (—E al ) as the canonically conjugate momenta of A“. In this 
way Eq. (1731) can be seen as a Hamiltonian density H(A, E) obtained from "Hnanm ($, n) via the substitution of the 
phase-space transformations 


(<fr,n )^(A,E), 


(82) 
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which follow from the inverses of Eqs. (l66l) and (l73l) plus Eq. (l68l) . 


4 | 

yj A b 3 A\ 



(83) 


in terms of the new variables. But, since the transformations are generated by the change of variables m 
in coordinate space, we know from classical mechanics that the full transformation in phase space is a canonical 
transformation. In this way we automatically recover the PBs 

{A“(x),^(y)} = 0, {£?“(x),^(y)} = 0, {A“(x),£^(y)} = -<5“ 6 ^<5 3 (x - y) (84) 


from Eq. m- To summarize, from each Hamiltonian version of the NANM, defined by the different values of n 2 , we 
can regain (via a canonical transformation) the Hamiltonian density d79l) together with the canonical algebra (1A29I) . 
The Hamiltonian density (1771) differs from the YM Hamiltonian density (IA28I) only in the fact that the Gauss laws 
fl b = ( DiE b — J h0 ) = 0 do not appear as first-class constraints, because Ag are not arbitrary Lagrange multipliers, 
but rather as functions of the coordinates, as shown in Eq. (l83l) . 


C. The evolution of the Gauss functions 


The calculation follows the same steps dictated by Eq. (1531) in the case of the SL-NANM. In the parametrization 
m, the result is 

= -gC abc A b 0 n c - + D 3 - D 3 (^4^) + A {^ A > b ) > ( 85 ) 

where Aq is given by Eq. (1831) . The above evolution equation leads to the same statements regarding the equivalence 
of the SO (N) YM theory with the NANM as those stated after Eq. (1571) in the case of the SL-NANM. 


IV. SUMMARY AND FINAL COMMENTS 

The possible interpretation of gauge particles (e.g., photons and gravitons) as the GB modes arising from some 
spontaneous symmetry breaking is an interesting hypothesis that would provide a dynamical setting for the gauge 
principle. 

In this paper we have taken the Nambu approach, whereby spontaneous SLSB is incorporated in an effective way 
in the model by means of a nonlinear constraint. These models can be understood as generalizations of the nonlinear 
sigma model describing pion interactions. The challenge posed by this setting is to show the conditions under which a 
the violations of Lorentz symmetry and of gauge invariance (introduced by the nonlinear constraint) are unobservable 
in such a way that the Goldstone bosons that appear can be interpreted as the gauge particles of an unbroken gauge 
theory. In other words, one tries to determine the conditions under which the corresponding Nambu model is equivalent 
to the unbroken gauge theory. Such conditions have been studied using perturbation theory for electrodynamics and 
YM theories, for example, in Refs. dam EMI. The main result is that, to the order considered (usually the 
tree-level or one-loop corrections) and after imposing the Gauss laws plus current conservation, the violations of 
Lorentz symmetry are unobservable, so that the corresponding Nambu model reproduces the corresponding gauge 
theory, with the gauge bosons realized as the corresponding Goldstone bosons. 

In this work we have generalized the nonperturbative Hamiltonian analysis developed for the Abelian Nambu model 
in Ref. [2^| to the non-Abelian case. Also, we have made an important conceptual and practical improvement in the 
method of dealing with the relation between the NANM and the corresponding YM theory. On the other hand, no 
discussion is provided here about the possible perturbative equivalence of the corrected Nambu model and the YM 
theory in a fixed gauge. 

In Sec. II we considered the specific case of the SL-NANM by solving the nonlinear constraint © in terms of A a l jzJ 3 
and starting with the remaining 4 N — 1 d.o.f. per point in coordinate space. The Lagrangian equations of motion were 
obtained, yielding different results than from the YM equations of motion, as expected. The canonical momenta and 
the canonical Hamiltonian were subsequently constructed, with the appearance of 2 (N — 1) second-class constraints. 
The standard variables of the YM theory, A“ and E b \ i = 1,2,3, were written in terms of the canonical variables of 
the SL-NANM, the canonical algebra of which induces the standard YM algebra for the former variables at the level 
of PBs. Appendix B includes a summary of the required PBs which prove the previous statement. The second-class 
constraints were further strongly imposed by introducing Dirac brackets, whose values for the variables A“ and E b i 
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turned out to be the same as the previously calculated PBs. The final extended Hamiltonian for the SL-NANM, 
rewritten in terms of the variables A“ and E b i , has the same form as the standard YM Hamiltonian, except that the 
Gauss laws fl a = 0 do not appear as first-class constraints. The time evolution of the functions fl a , according to the 
SL-NANM dynamics, were calculated, yielding the result that after demanding current conservation, the imposition 
of the Gauss laws at some initial time yields fl a = 0 for all times. The final 3 N d.o.f. in coordinate space of the 
NANM were recovered, since \ (2(4 N — 1) — 2 (N — 1)) = 3 N. It was emphasized that a similar analysis could be 
carried out for each of the remaining cases corresponding to the time-like and light-like versions of the NANM. 

Section III presented a substantial conceptual and practical improvement over the previous individual calculations 
to study the relation between NANM and YM theories. Starting from an alternative parametrization that solves the 
constraint © for arbitrary values of n 2 in terms of 3 N d.o.f. A = 1,2,3, we showed that the calculation of 
the NANM canonical momenta can be written in such a way that the chosen parametrization induces a direct 

relation between the YM variables A“ and E b i and the canonical variables of the NANM. Since the phase-space 
transformations are induced by a coordinate transformation (the chosen parametrization), we know from classical 
mechanics that the phase-space transformation ($( 4 , n^) — > (A“, E b ' J ) is a canonical transformation, once the E hj are 
recognized as the momenta canonically conjugate to the Af via the kinetic part of the NANM Hamiltonian action. In 
this way, one immediately concludes that the resulting algebra of the (A“, E b i) has to be the canonical one, without 
requiring the detailed and tedious calculations that were necessary in the discussion of the previous section. It is 
interesting to observe that the identification of the canonical transformation is independent of the detailed structure 
of the chosen parametrization, as soon as it provides an invertible change of coordinates $(4 = $(4 (A“). Another 
useful property of this parametrization is that it exhibits the NANM as a regular theory (i.e., no constraints appear in 
the Hamiltonian analysis), since it includes just the necessary 37V d.o.f. of the NANM. This was proved in Appendix 
C. The canonical NANM Hamiltonian, rewritten in terms of the YM variables (A“, A 6 - 7 ), again has the same form as 
the YM Hamiltonian, except that the Gauss laws do not arise as first-class constraints. The time evolution of the 
functions was also calculated, with similar results as in the SL-NANM. 

The relation between our approach and the method of Ref. []j|, which also included pure Goldstone field modes, 
was elucidated in the paragraphs after Eqs. (l55l) and (157>1) . 

To summarize, a nonperturbative equivalence between the SO(IV) YM theory and the corresponding NANM has 
been established, after current conservation and the Gauss laws are imposed as initial conditions for the latter. 
Actually, the Gauss laws — valid now for all times — are next added as Hamiltonian constraints —N a tt a to Eq. (1791) . 
with arbitrary functions 7V°. The further redefinition Ag + N a = 0“ leads to the final YM Hamiltonian 

H ym = ^(E 2 + B 2 )-0 a n a + J“A“, (86) 

where 0 a are now arbitrary functions. In other words, the Gauss laws are imposed a la Dirac upon the physical states 
IH 7 )phys by demanding that H a |4') p hy S = 0. Also, the emergence of the SO (N) YM theory subsequently guarantees 
current conservation for all times, as a consequence of the YM equations of motion. The established equivalence is 
independent of any gauge fixing and supports the idea that gauge particles arise as the Goldstone bosons of a model 
exhibiting a spontaneous Lorentz symmetry breaking that is not physically observable. 
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Appendix A: The SO(N) Yang—Mills theory 


We present a brief review of the Hamiltonian formulation of the standard SO (TV) YM theory. The main motivation, 
besides establishing some notation, is to recall the basic properties of the YM theory that have to be recovered in 
order to state its emergence from the different versions of the NANM. 

The Yang-Mills Lagrangian density is given by 


C = Tr 



_ J 7 ‘A 


M ) 


(Al) 


where boldfaced quantities denote matrices in the Lie algebra of the internal symmetry group SO (TV) with N(N — 1) /2 
generators f a ; i.e., M = M a t a . This algebra is generated by [t a ,f b ] = C abc t c , where the structure constants C abc are 
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completely antisymmetric. The field strength is 


E 

A fl V 

— T g [A m , Ajy], 

(A2) 

and the equations of motion are 

d^ v = r, 

(A3) 

where the covariant derivative is defined as 

D„M = d,M +g [A m ,M], 

(A4) 

From the above definitions we obtain 

[Dp,D v ] M= [M,F„„], 

(A5) 


which leads to current conservation, D„3 1 ' = 0. 

The expressions of Eqs. (IA2I) and (1X41) in terms of the components of the corresponding fields are 

F; v = d»K - d v Al + gC^A^Al, 

(D^M) a = D^M a = d^M a + gC abc A b ^M c . 

The Jacobi identity for the connection is 

(jarb£jcdr _|_ ^jcrb^jdar _|_ ^jdrb^jacr q 


(A6) 

(A7) 

(A8) 


in terms of the structure constants. The group indices a = 1,2,..., N are raised (lowered) by the metric 6 ab ( 6 a b ) and 
their position as superscripts or subscripts is just a matter of convenience in writing the corresponding expression. 
Next we review the Hamiltonian version of the YM theory. The canonical momenta are given by 


H°m _ 


dL 

aW 


Therefore, considering 


we find 


d(F^Fr) 

d(A%) 


4p0a 

a > 


ng = o, 


H“ = F? 0 = -E9, 


which satisfy the nonzero PBs 

(Ag(x,t),n^(y,t)} = 6 ab S(x — y), {A?(x,t),n bj (y,t)} = <5^ ob J(x- y). 


(A9) 

(A10) 

(All) 

(A12) 


In the following we assume that all PBs are calculated at equal times and we suppress the label t in most cases. From 
Eq. (IA6I) we get Af as 


Af = E? + d.A* + gC abc A b A c 0 = E“ + D 


(A13) 


We also introduce 


B a k = 2^ kF%. 


(A14) 


Recalling that 


D^(N a M a ) = (D^N a ) M a + N a (D /1 M a ) = d fJ ,(N a M a ), a = 1,2,..., N, (A15) 

which allows us to perform integration by parts within the action, we find the canonical Hamiltonian density 

1 


H = K A t - £ = 7j( E + B ) - A o(F>iEi - -n a - Ji A “, 


(A16) 
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where 

E 2 = tr( F 0l F 0l ), B 2 = ±ir(F y F«). (A17) 

We employ Dirac’s method to construct the canonical theory, due to the fact that primary constraints 

E a = ng ~ 0, (A18) 

are present. The extended Hamiltonian density is given by 

Ue = \{ E 2 + B 2 ) - A&DiEi - J°) a + J^A\ + A a E a , (A19) 

where X a are arbitrary functions. The evolution condition of the primary constraints 

E“(x) = {E°(x), J d 3 y H E ( y)} ~ 0, (A20) 


leads to the Gauss laws 


fl a = (DiEi - J°) a ~ 0. (A21) 

It is not difficult to prove that Eqs. (IA18I) and (IA21I) are the only constraints present and that they constitute a 
first-class set. In fact, calculating the time evolution of f l a yields 

f2 Q (x) = j^ Q (x), J d 3 y H E (y )} = {^“(x), J d 3 y Qb 2 - A b Q b + J9A^j } - d 0 J 0 “■ (A22) 

From the PBs (IA12I) together with Eq. (IA 8 I) we obtain 

Jd 3 y{n a (x), B 2 (y)}=0. (A23) 

The PB of the constraints (1A21I) produces 

|A £?(x), J d 3 y D 3 £‘(y)J = C abc D k E c k (x), (A24) 

which leads to 

J d 3 y (B a (x), B b (y)} M b ( y) = C abc M b [B c + J 0 C ] (x). (A25) 

In this way, 

n a = -do JS - C abc A b JS - C abc A b n c - D k J ak = -C abc A b n c - (A26) 

which is zero, modulo the constraints and using current conservation. 

Normally one fixes 

ng ~ 0, Ag ~ 0 a , (A27) 

with 0 a being arbitrary functions to be consistently determined after the remaining first-class constraints f l a are 
fixed. 

The final Hamiltonian density is 

He = i(E 2 + B 2 ) - 0 Q (AA - J°)a + J-K- (A28) 

Once ng and Ag are strongly fixed, the Dirac brackets of the remaining variables are 

{A“(x,i),A*(y,t)}* =0, {E ai (x,t),E bj (y,t)}* = 0, {A“(x, t), E b \y, t)}* = -6{5 ab S( X - y). (A29) 

The final count of the number of d.o.f. per point in coordinate space yields 

#d.o.f. = i (2 x AN - 2 x 2N) = 2 N. 


(A30) 
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Appendix B: The bracket algebra in the SL-NANM 


In this case the canonical variables are 
1, 2,3 and t = 1,2. The nonzero PBs are 


Ag, Ai, Af, big, n lfc , n ai where a = 1,2,...,TV; a = 2,3,...,7V; i = 


{^E5( x i*)> n o(y,i)} = <5 3 (x - y)S ab , {Ai(x, t), II lfc (y, t)} = S k <5 3 (x — y). 

{Af(x, t)} = <f?W(x - y). (Bl) 

The YM canonical variables (Af, E b ) are given by 


A\ = ^P + (Af ) f-(A\f-(A^ 2 , A\, Af, 


pi _ /ii | n o 
A, - A, I ^ 


Ef = n 11 + [ ) A, 


pa = jja 


55', a- 


(B2) 


in terms of the canonical variables of the SL-NANM. 

In the following, we do not specify the coordinate dependence of each term: it is to be understood according to the 
following convention 


{P(x),Q(y)} = {P,Q}. (B3) 

Additionally, we also suppress the unit 5 3 (x — y) in coordinate space. We do not provide any details for each 
derivation: we only include the final results. 

Our first goal is to calculate the equal-times algebra among the YM variables Af, E b defined in Eq. (IB2I) in terms 
of the canonical algebra (IB1I) of the SL-NANM. To this end we first consider the PBs between the A 3 [which is the 
solution of the constraint (p[j)] and the canonical conjugate momenta of the SL-NANM. The results are 


{n Bfc ,Ai} = 


3 


{ng,Ai} = - 


A a 

—1 

3 


{n 


lfc 


1 M 


(B4) 


1. The A — A sector 


The algebra among the canonical coordinates (Ai, Af) and A 3 is trivial because the former satisfy the canonical 
relations and A 3 = Ag(Ai, Af), so that we have {Af, A b } = 0. 


2. The A — E sector 


We obtain 



{Al,El}=l, 

{Al,El}=0, {A\,E?}= 0, {A\,E \} = 0, 

(B5) 


{Al,E\}=5l {A\,El} 

= 0, {Af,E\} = 0, {Af,E\} = 0, {4,Pf} = S^Sf. 

(B 6 ) 

The above relations can be summarized as 





{Af, E k } = S ab 6l 

(B7) 



3. The E — E sector 


We have 

{El,Ef}=0, { E.{, Ef } 

= 0, { E l E l)=^ {El,El}= 0, [El,E b ) = 0. 

(B 8 ) 


To summarize, the PB algebra of the canonical variables in the YM theory, calculated from the canonical algebra of 
the corresponding SL-NANM, reproduces the YM algebra (IA29I) . 
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4. The (AS, ns) - (A?, *£) sector 


We also need the PBs between the Ag, IIg and the YM variables Ag, The results are 


{Ag, A b }=0, 


K, E b } = 


(B9) 


{ns, A b } = -a 16 a?^| 
{ng,^} = §^ 


*3 0 


{ng, £{} = § 


4« 41 

jal _ ^0^0 


(^r 


{ns, 4} 


^Aj d ' 


(BIO) 


5. The - (Ag, Ag, £“) sector 

The equality of the Dirac brackets and the PBs derived in Eq. (15U1) is a direct consequence of the PBs between the 
constraints = Ilg — IljAg/Aj and the YM variables, which are 

{n,E\}= 0, {$?,£} }=0, {<S>l,E\ } = 0, }=0, 

{$“, AJ} = ^, {$?, A* } = -,5 s5 , {$?, A 3 } = 0 , {$?, A)} = 0 , 

{$?, A b }=0, {$?, A£}=0, {$“, Bj}=0. (Bll) 


6. The constraints sector 

Next we provide the results for the calculation of the PBs of the constraints and 

{$“, $51 = Q 


The calculation of 


<J>^j requires 

{$“, n b }=o, 

in virtue of the relation {<!>“,£{}= 0 calculated in (IB11I) . The final result is 




where we also have made use of the constraint fUAg = AgfU from Eq. 


Using the results 


we calculate 


(Ag, D b }=0, {Ag, D b }=0, 


{$g, $ b } = (AJ) 


n 2 J n a n b 
u 1 ’u 1 


The last PB is 


f D a 

n b ) 

[ 1 ( 


n 1 j 

m 2 \ 


r>b na 

s~iabm 1 s~ilam s^ilbm 

6 +c m~ c m 




(B12) 

(B13) 

(B14) 

(B15) 

(B16) 

(B17) 
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where the relation 


{Q a , Sl b } = C abc DiEZ, 


(B18) 


has been used. The final result is 

$ 2 } 


( 3 . 

yn 1 2 3 


2 



Q b 

6 m “ c 



(DiE ™) • 


(B19) 


Appendix C: The Wronskian 

In the formulation of the NANM presented in Sec. Ill, the Wronskian arising from the Lagrangian is 

Wnanm = det 


f)A c f)A c ' 

* f), i,A = 1,2,3, a = l,..., IV, 


d& b B d$ a A 


where the coordinate transformation Af = Af{<& b B ) is invertible so that 


det 


dAf 

d 5^ 


^0. 


(Cl) 


(C2) 


We next show that the property 


coordinate transformation from 3N Af to 3IV $ b B as 


guarantees that Wnanm is nonzero. To this end, it is simpler to relabel the 


Ar = A h (<3>s), det 


\d$s) 


7^0, R,S = 1,2,3,. ..31V 


The required Wronskian is 


, JdA R dA R \ 

Wnanm = det ■ 


In terms of the 3IV-dimensional epsilon symbol e T ' T2 ’"' ,T:iN , the above can be written as 

rn _ 7'i I j..... 7:;.\ ( dA Rl 3 Ar 3 ^ / 0Ar 2 dA R2 ^ / BAr 3N 0Ar 3N 

NANM “ 6 IWl d® Tl J V ) "■ V d®3N a$T 3JV 

_ f SAr 1 8Ar 2 8Ar 3N \ { Ti t 2 ,...,t 3N dA R 1 \ 


Wnanm 


y 3 $ 2 d$>3N 


d&n 9 ®t 2 d®T 3 N J 


But 


so that 


( f TiTa . T 3N dA Ri dA R 2 0 Ar 3N \ _ Rl R 2t ... : R 3N , , ( 9Ar \ 
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det 


\d<S> T )_ 


± 0 , 


employing 
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